A PHYSICAL EXPLANATION FOR THE TILDE SYSTEM IN 

THERMO FIELD DYNAMICS 

DONG MP, HE-SHAN SONG^ 
Department of Physics, Dalian University of Technology, Dalian 116024, 

P.R.China 

^ . YING AN 

O . 

, Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 

H . 

; 100080, P. R. China 

For a two-body quantum system, any pure state can be represented by 

m . 

CSJ , a biorthogonal expression by means of Schmidt decomposition. Using this 



in the composite system which include a thermodynamic system and its sur- 



(N 
> 

, roundings, it is found that the tilde system in thermo field dynamics is just 



^4 , the surroundings of the real system. 



PACS numbers: ll.lO.Wx, 05.30.-d, 03.65.Yz 

In 1975, Y. Takahashi and H. Umezawa constructed a quantum field theory at finite 



C3 

\ temperature: Thermo Field Dynamics (TFD)^. In this theory, the "fictitious" (or tilde) 
^3 ! 

system is essential to convert thermal statistical average into expectation value with respect 
^ ^ to a pure state. The central idea of TFD is the doubling of the Hilbert space of states. 
I The operators on this doubled space are effectively doubled as well^'^. After that, most of 
the techniques in quantum field theory can be effectively used to deal with the problems 
in statistical physics. However, so far there exists a question in TFD which has not been 
answered satisfactorily: what is the physical meaning of the tilde system? Is it really a 
fictitious system? The discussion of the physical meaning for the tilde system started from 
the birth of TFD. For example, Takahashi and Umezawa suggested in their original paper 
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that the particles in the tilde system may be regarded as the holes^. Y. X. Gui has suggested 
it may be regarded as the mirror world of our real one^. Before the passing away of Umezawa, 
he still referred to the tilde system and tried to explain its physical meaning. He pointed 
out that the tilde operator may be called the thermal operator for all the thermal quantities 
need the tilde operators^. We think these explanations are ambiguous and should belong to 
guesswork without good ground. In this letter, we first consider a two-body quantum system, 
its any pure state can be represented by a biorthogonal expression by means of Schmidt 
decomposition. Notice that any thermodynamic system is not a truly isolated system, it is 
always in contact with the surroundings. Using the idea about Schmidt decomposition in 
the composite isolated system which include a thermodynamic system and its surroundings, 
we find that the tilde system in thermo field dynamics is just the surroundings of the real 
system. 

Firstly, let us imagine an isolated system with an definite state vector] ■0). This isolated 
system will be supposed a composite system and to consist of two parts. One is the sys- 
tem whose properties we wish to study. The remainder, to be called the surroundings of 
the system, will not be observed. The system to be measured has a complete set of com- 
muting observables whose eigenvalues will be labeled collectively as a Roman letter, as i. 
The surroundings have a complete set of commuting observables whose eigenvalues will be 
labeled collectively as a Greek letter, as /i. The observables belonging to the system and 
surroundings commute because they refer to different things, so a basis for the isolated sys- 
tem, system (denoted by A) plus surroundings (denoted by B), is the set of simultaneous 
eigenstates^ \i) <^ |//) = For simplicity, it will be assumed that all the eigenvalues are 

discrete, so the completeness relation is 



(1) 



By using equation (1), j'^) can be expanded as 



(2) 
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with the expansion coefficients Oj,^ = . 

An observable Fa which refers to system A alone, when considered as an observable of 
the composite system, should have the following form 

F^Fa® Tb, (3) 

where Ib is identity operator which acts on the Hilbert space of the surroundings. The 
expectation value of Fa <S) Ib in the state 1-0) is as usual 



Fa ® Ib) = {i^ Fa^ Ib ^) ■ (4) 
With equation (2), we can write this as 

'Pa ®Ib) =Y1 '^*3,u ^\ \Fa ®Ib\Y^ \h /«) 

= tTA (paFa) , (5) 



where 



= tvB (pab) (6) 

is called the reduced density operator (or matrix) of the system^. 
It is easy to see that pA has the following properties 

PAt = Pa, (7) 
tVAPA =1, (8) 

Pa is semi-positive definite. (9) 

Since the density matrix pA is a semi-positive definite Hermitian matrix, it can be diagonal- 
ized and it will be have non-negative real eigenvalues. Notice that the relation p\ — pA, in 



general, does not hold. This means that the system is in a mixed state. Let the set {\i >} 
be the complete orthonormal eigenstates of the density matrix pA and {Pi} its eigenvalues 
{Pi > since pA is semi-positive definite), then the density matrix can be written in the 
form 

p^ = ^p,K)(z|. (10) 

i 

where Pj is the probability of finding the system in state \i) and Z^Pj = 1. In this way, the 

i 

pure state 1-0) of the composite system can be rewritten as the following form 

m^EPi^'\h^^)^ (11) 

i 

where 

i«.) = E^r'^'«^,Mi/«)- (12) 

Here, {!«,>} is a set of complete orthonormal states of the surroundings. In fact, equation 
(11) is the Schmidt decomposition of a pure state of the composite system^, and the set 
ccj >} is called Schmidt basis. Prom the viewpoint of the surroundings B, a pure state 
of the composite system (A+B) appears as a mixed state of the system A, described by 
a density matrix obtained by tracing over the degrees of freedom of the surroundings, and 
the density matrix pA is diagonal in the Schmidt basis. That is equation (10) can also be 
obtained by the following way 

Pa = trs (Pab) ^ trs {ip\) 

= E^^K)(^I- (13) 

i 

We denote the expectation value of a quantity in mixed state by the symbol (), for this 
expectation value is a double averages including the thermo and quantum averages. Thus 
we can express equation (5) as 

{Fa^Tb) = {P^. (14) 
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Using equation (11) and (10) or (13), the above equation can be written in the Schmidt 
basis as 

i,j i 

Therefore, we can see that the expectation value of a variable Fa of the subsystem A in mixed 
state can be converted into the expectation value of a variable Fa ® Ib of the composite 
system with respect to a pure state. As a result, it is inevitable to enlarge the original 
Hilbert space of the system to include the Hilbert space of the surroundings, and hence the 
corresponding degree of freedom is doubled. 

Now, consider a thermodynamic system in thermal contact with a heat reservoir. We 
know that a macroscopic system to be dealt with in statistical physics consists of enormous 
number of microscopic elements interacting in a complicated fashion with the surroundings, 
the energy intervals between adjacent eigenstates will be much smaller than any distur- 
bance due to the surroundings. In other words, any small disturbance will cause transitions 
between eigenstates, resulting in a uniform probabilistic distribution. Consequently, the be- 
havior of the dynamic system may be expected to be described by certain probability laws^. 
This implies that the state of a system in equilibrium may be regarded as an incoherent su- 
perposition of eigenstates of the system. As is well known, a collection of identical systems 
used for studying probability characteristics is called an ensemble. Therefore, the states of a 
thermodynamic system that are described by statistical ensemble are mixed states, which are 
distinct from the pure states which are described by state vectors in quantum mechanics^. 
In short, owing to the uncertain interaction between the thermodynamic system and its 
surroundings, an uncertainty about micro-state is produced. In the way, we know only an 
ensemble of possible states, say \n >, and the probabilities p„ of the system being in state 
|n >, respectively. The information on the system is given by the density matrix. 

In statistical physics we always deal with systems that interact with the surroundings. 
Here we can regard the system plus its surroundings as a truly isolated system. Following 
the former mark, we denote the thermodynamic system and its surroundings by A and B, 
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respectively. Suppose Fa is an observable of the thermodynamic system. As is well known 
that the thermal statistical average of this observable in the thermal equilibrium state is 
given by 

{pAj^Tr{pAFA), (16) 

where, pA is the density operator of the thermodynamic system; Tr -Fa) denotes the 
trace of the operator PaFa and is the sum of all the diagonal matrix elements of in any 
representation. 

For definiteness, consider now a closed thermodynamic system with the numbers of 
particles, volume and temperature fixed (The following discussion can be easily generalized 
to the open system which may be described by grand canonical ensemble.). Let Ha be 
the Hamiltonian operator of the system, then the density operator of the system under 
consideration has the following form^'^ 

Pa = (17) 

where 

Za = Tr (e-^-^^) (18) 

is the partition function of the system, and /3 = (kBT)"^ , with is Boltzmann constant, 
T the temperature of the system. 

Notice that the density operator and the Hamiltonian operator of the system commute, 
they possess a complete orthonormal set of simultaneous eigenstates. So, both pA and Ha 
are diagonal in energy representation. Assume we have the following eigenvalue equation of 
the Hamiltonian 

HA\n)=E^\n). (19) 

A set of linearly independent eigenvectors \n) of Ha can be chosen to be orthogonal and 
normalized 
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{m\n) = 5mn, (20) 
and the completeness relation can be express as 

Y,\n){n\^lA. (21) 

n 

Then, the equation (17) becomes 

pA^^Pn\n) {n\, (22) 

n 

where 

is the probability for the system to have the energy eigenvalue Em and it satisfies the 
normalization condition X^Pn = 1- Thus, the thermal statistical average of the quantity Fa 

n 

can be expressed as 



FA) = J2Pn{n\FA\n). (24) 

n 

Suppose the statistical average of a quantity Fa can be written as expectation value 
with respect to a pure state |0(/3)). According to equation (14), we can write the following 
relation 



'Fa) = {0{P) Fa®Ib 0{(3)). (25) 

According to equations (2), the pure state |0(/3)) can be expanded as 

|0(/3)) = ^a„,Jn,/x), (26) 

where a„^^ = ('^^A*! |0(/3)), and the sets {|n >} and {|// >} are the complete orthonormal 
eigenstates of the thermodynamic system and its surroundings, respectively. 

In the energy representation, by means of Schmidt decomposition, \0{(3)) can be written 

as 
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\Om^Y.pT\n.o.n). (27) 

n 

with 

Obviously, the equation (25) automatically hold, 

X] Pn {n\ Fa \n) = ^ pH^pH'^ {m, q;^| Fa ® Ib \n, a^) . (29) 

n m,n 

Now, we denote q;„ by n, the pure state |0 (/?)) can be written as 

\Om-Y.pT\n.n) (30) 

n 

where p„ is given by equation (23). Equation (30) is nothing but the thermal vacuum state 
which was introduced in thermo field dynamic by Takahashi and Umezawa^. So, we can see 
clearly that the tilde system in TFD is just the surroundings of the real system. 

In addition, we can see from equation (30) that the thermal vacuum state \0 {(5)) is an 
infinite superposition of the direct product state \n) ® \n). In general, the pure state \0 {[3)) 
of the compound system cannot be expressed as a product \n)®\n) of pure states of its part 
(i. e., A and B) and hence is an entangled state^^. We notice that, for the system under 
consideration, the pure state \0 {(3)) is a partly entangled pure states, for the coefficients 

1 /2 

{pn) are not equal to one another. Obviously, this conclusion is true for an open system 
which may be described by grand canonical ensemble. In addition, for an isolated system in 
statistical physics (In fact, it is not a truly isolated system.), p„ in equation (30) becomes 
a constant. Since all of the expanding coefficients in equation (30) are the same now, the 
state |0 (/?)) corresponding to an isolated system is a maximally entangled state-*^^. 

The density matrix of the compound system {A+B) corresponding the state vector 
\0{I3)) is given by 

PAB = \0 m^BAB {O (/?)| = E pU'pI/" H ® (^1 ® (^1 • (31) 

m,n 
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Because the |0 is a pure state of the compound system (A+B), the density matrix pab 
is non-diagonal. However, we are only interested in the real system A. When we measure 
a dynamical variable of the real system, say energy, the state of the system will transfer to 
one of its energy eigenstate and lead to the quantum decoherence^^. Mathematically, this 
means the density matrix will be traced in the tilde states. That is we should describe the 
real system by means of the following reduced density matrix^'^ 

PA = Trs (Pab) = E P^b |^ = E PrnV^' {l\ \m) |m) {n\ {n\ ||) = \n) {n\ , 

I l,m,n f- 

(32) 

which is the same as the corresponding one in statistical physics (see equation (22) ). So, the 
surroundings forms an pure entangled states with the thermodynamic system upon tracing 
out the former we can get the correct mixed state for the latter system. 

In summary, by the unavoidable interaction of the thermodynamic system in question 
with the surroundings, the thermodynamic system will be in a mixed state. If the thermal 
statistical average is replaced by a vacuum expectation value in a pure state, then the 
Hilbert space of the thermodynamic system should be generalized to include the one of the 
surroundings, and the corresponding operator become a direct product which include the 
operator of the system and the identity operator of the surroundings. 
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